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Abstract. We present the results from global fits of inclusive DIS experimental data using 
the Balitsky-Kovchegov equation with running coupling. 



On rather general physical grounds, one expects that the large gluon density that develops 
in the proton (nucleus) wave function with decreasing Bjorken-x will lead to deviations from 
standard collinear factorization. The necessary unitarity of the theory requires the rate of growth 
of the gluon density to be limited. Such a taming arises naturally once non-linearities, due to 
gluon overlap, are accounted for in the evolution. 

The Colour Glass Condensate (CGC) 1-3 1, along with its B-JIMWLK evolution equations 
[4j[5] , provides a solid theoretical framework in which to address the non- linear dynamics driving 
the small- x evolution of the proton (nucleus) wave function. The mathematical complexity of 
the CGC evolution equations has hindered its direct phenomenological application. In face of 
such difficulties, the phenomenological use of the CGC has proceeded via dipole-models in which 
the general physical properties of the CGC have been, in one form or another, implemented in 
an effective manner. The behaviour at small- £ of the parton distribution functions plays an 
important role in the computation of many observables in proton-proton, nucleus- nucleus, and 
proton-nucleus at the LHC. Given that evidence for small, yet systematic, deviations from NLO- 
DGLAP evolution have been identified in HERA DIS data [6], it is mandatory that the best 
available and numerically tractable non-linear small- x evolution tools replace CGC-motivated 
models. 

The simplest theoretically sound realization of the CGC, the Balitsky-Kovchegov (BK) 
equation (5} [7], although tractable from the numerical point of view, fails to reproduce 
experimental data unless rather unphysical assumptions are made. Early heuristic modifications 
of the BK equation |8| indicated that NLO corrections would drive the evolution towards 
compatibility with experimental data. The explicit computation of the full NLO-BK [9]-|ll| 
paved the way to the rigourous and phenomenologically relevant use the CGC. Although these 
corrections present a rather complicated structure and are not presently amenable to an efficient 



numerical implementation, it was shown [12] that considering only running coupling corrections 
to the BK equation accounts for most of the NLO effects. The ability of the rcBK equation 



to correctly describe experimental data was shown in 13 . There, a global fit of the inclusive 
DIS structure function led to the public release of a parametrization of the dipole scattering 
cross section at small values of x. Thus, running coupling BK (rcBK), the most reliable and 
phenomenologically usable non-linear small- £ evolution tool, has become the standard choice 



for CGC computation of experimental observables (e.g., 14-16]). 

The HI and ZEUS combined data set for the DIS reduced cross section [T7] poses a series 
of new challenges and opportunities. On the one hand, the high-accuracy of the data puts 
the rcBK equation under much more stringent test conditions and will eventually lead to a 
better constrained parametrization. On the other hand, the data being given directly in terms 
of the reduced cross section, the experimentally measured quantity, eliminates the theoretical 
bias associated with the extraction of F 2 and Fl from data. Here, we present the first results 
obtained both from the global fit of these data within the setup of [13] where only light quarks 
were taken as contributing to the DIS cross section, and from a global fit in which available F% 
data is used to constrain the heavy quark contribution to the DIS cross section. 

In the kinematical domain (Q 2 < 50 GeV, x < 0.01) considered in this work, the reduced DIS 
cross-section a r can be written in terms of the virtual photon- proton cross-section gt l, with T 
(transverse) and L (longitudinal) the polarization of the virtual photon, as 

( tf\ Q 2 I , 2(1 -y) \ 
^ VM) = 4^a— [ aT + l + (l-yy aL ) ■ (1) 

Here, y = Q 2 /sx (-^/s is the center of mass collision energy) is the inelasticity variable. In the 
dipole formulation of QCD, valid for small-x, one has 
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where ls ^ ne light-cone wave function for a virtual photon to fluctuate into a quark- 

antiquark dipole of quark flavor / (with mass rrif and electric charge ej). M(r,x) is the 
imaginary part of the dipole-target scattering amplitude averaged over impact parameter, with r 
the transverse dipole size, aoj is (half) the transverse area over which quarks of a given flavour 
are distributed. Light quarks are taken to be identically distributed (<7o,/=«,<i,s = o"o,/ig^)- When 
accounting for heavy flavour contributions in ([!]) we allow o"o,/= Ci 6 = &o,heavy to be different from 
&o,Ught- These are free fit parameters. 

The evolution of A/*(r, x) is given by the rcBK equation: 
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where r2 = r — ri, and xq (in our case xq = 0.01) is the value of x where the evolution starts (the 
highest value of x in the data included in the fits). The kinematic cost of the virtual photon to 
fluctuate into a dipole of a given mass is accounted for by x = x exp (1 + 4m 2 /Q 2 ). The running 
coupling in ^ is evaluated at 1-loop accuracy in coordinate space 



where the constant C (a fit parameter) accounts for the uncertainty in the Fourier transform from 
momentum to coordinate space. When only light flavours (/ = u, d, s) are taken as contributing 
to ([!]), the coupling ([4]) is evaluated with rif = 3 (light only). Once the contribution of heavy 
flavours (/ = c, b) is accounted for, the coupling Q is evaluated in a variable flavour scheme 
matched on the quark mass thresholds (light+heavy) . Since in the rcBK equation ^ all dipoles 
sizes are explored, an infrared regulation is called for: a s (r 2 > r 2 r ) = aj r = 0.7 where r? r is 
the dipole size at which the coupling reaches aj r . Although we have explored |18| a number 
of different initial conditions, the results show a negligible dependence on specific choice. For 
brevity, only results obtained with the GBW initial condition 



M GBW (r,x = x ) = 1 - exp 
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(5) 



are shown here. The fit parameters Q 2 and 7 are, respectively, the saturation scale and the 
characteristic fall-off of the dipole scattering amplitude with decreasing r at the initial xq. 

Fig. [l] shows a comparison of data for the reduced cross section (selected Q 2 data bins) with 
the results of a light only fit using all combined Hl/ZEUS data (196 points) within cuts Q 2 < 50 
GeV 2 and x < 0.01. The mass of the light flavours was kept fixed, mu g ht = 0.14 GeV. The 
obtained fit parameters — (JQ,u g ht = 33.1 mb, Qg tag/it = 0-^6 GeV 2 , "yn g ht = 0.97, and C = 2.3 
- give x 2 / d -°- f - = 206/192 = 1.07. 
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Figure 1. Comparison of 
reduced cross section experi- 
mental data (in black) with 
our (light only) fit results (red 
stars). Notice that the exper- 
imental error is often smaller 
that size of the symbols used 
for the data points. 
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Fig. [2] shows a comparison of data on (the charm contribution to F2) with the results 
of a light+heavy fit. In this case, all existing data on F£ (within the same cuts as above) is 
added to the Hl/ZEUS sample (a total of 265 points). We allow for independent light and heavy 
initial conditions and for a free light flavour mass. The heavy quark masses are set to m c = 1.4 
GeV and = 4.75 GeV. The x 2 /d.o.f. = 312/257 = 1.2 is obtained for the fit parameters 
CT ,u ght = 32.7 mb, Q 2 ^ Hght = 0.26 GeV 2 , llight = 1.24, a 0Mavy = 28.4 mb, Q 2 s0Mavy = 0.18 
GeV 2 , lheavy = 0.95, C = 3.3, and m light = 0.104 GeV. 

The ability of the rcBK equation to accurately describe the high precision HI /ZEUS combined 
sample and the charm component of the structure function strengthens its role as a viable 
phenomenological tool to address the small-x behaviour of the proton structure. Further details, 
along with a parametrization of the dipole-proton cross section for use in the computation of 



physical observables, will be given in 18 . 
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Figure 2. Comparison of 
i*2 data (in black) with our 
(light+heavy) fit results (red 
stars), light+heavy fits give 
description of reduced cross 
section similar to that shown 
in Fig. [T] 
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